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Abstract 

There are a number of algebraic classifications of spacetimes in higher 
dimensions utilizing alignment theory, bivectors and discriminants. Pre- 
vious work gave a set of necessary conditions in terms of discriminants 
for a spacetime to be of a particular algebraic type. We demonstrate the 
discriminant approach by applying the techniques to the Sorkin-Gross- 
Perry soliton, the supersymmetric and doubly-spinning black rings and 
some other higher dimensional spacetimes. We show that even in the case 
of some very complicated metrics it is possible to compute the relevant 
discriminants and extract useful information from them. 



1 



2 



A. Coley, S. Hervik, M. Durkee and M. Godazgar 



1 Introduction 

Recently there has been considerable interest in the study of general relativ- 
ity (GR) in higher dimensions and, in particular, in higher dimensional black 
holes [T]. The underlying motivation for this comes from supergravity, string 
theory and the gauge-gravity correspondence, but higher-dimensional GR has 
developed into a field of interest in its own right. 

It has been shown that even at the classical level, gravity in higher dimen- 
sions exhibits much richer structure than in four-dimensions (4D). One of most 
remarkable features of 4D GR is the uniqueness of the Kerr black hole, but this 
fails in a striking way in higher dimensions. There now exist a number of dif- 
ferent asymptotically flat, higher-dimensional vacuum black hole solutions [1], 
including Myers-Perry black holes [2] , black rings [S] |4] , and various solutions 
with multiple horizons (e.g. [5J|H]). 

The algebraic classification of spacetimes has played a crucial role in under- 
standing black holes in 4D (for example, in the discovery of the Kerr metric, and 
the study of its linearized perturbations U]). Algebraic classification in 4D can 
be described in several different ways, using null vectors, 2-spinors, bivectors or 
scalar invariants. There are strong links between these approaches; indeed, each 
of them can be used to give a different description of the 4D Petrov classification 
scheme. Algebraic classification has been generalized to higher dimensions using 
each of these different methods [8j [9j [10] , but it turns out that each approach 
leads to a distinct classification in higher dimensions [TTl [TU] . 

For vacuum solutions, where the Ricci curvature tensor vanishes, algebraic 
classification reduces to the classification of the properties of the Weyl curvature 
tensor. The most comprehensive and (to date) well-studied approach [^ [T^ITB] 
classifies null frame components of the Weyl tensor according to their boost 
weights under local Lorentz boosts by identifying a choice of null frame such 
that components of high boost weight vanish. More explicitly, one works in a 
particular frame 

{£,n,m,}, i^2,3,...,D-l (1) 

where £ and n are linearly independent null vectors, transforming as 

£ 1-^ X£, n 1-^ X~^n, nii i— rrii (2) 

We will refer to this as the alignment classification. This classification divides 
spacetimes into 6 primary different types: G, I, II, III, N, O (for definitions and 
further details of these various algebraic types see, for example, [51 1131 HH US] ) ■ 

Examples of spacetimes that are algebraically special in this classification 
include the type D Myers- Perry black holes (i.e., they admit two independent 
null vector fields with the type II property). Black ring spacetimes are also 
algebraically special on the horizon [1^; this will be discussed in more detail 
below. Applications of this classification scheme include the recent construction 
of gauge invariant variables describing perturbations of any spacetime that is of 
type II or more special [171. 

Compared to four dimensions, the algebraic types defined by the higher- 
dimensional alignment classification are rather broad, and it has proven more 
difficult to derive general results. Therefore, it is natural to ask whether there 
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is a useful way of refining these types further. The higher-dimensional bivector 
classification [TU] achieves this by analysing the bivector map 

C : Xfj^^ i~> ^j^''" X per , (3) 

which maps the space of spacetime bivectors (2-forms) X to itself. In four 
dimensions, a decomposition into (anti-)self dual parts can be used to render the 
bivector classification identical to the standard Petrov algebraic classification]^. 
However, in higher dimensions this is not possible. 

In TU^ , two of the present authors analysed the bivector map in higher di- 
mensions, and showed how its eigenvalue structure is related to the alignment 
classification [9]. In particular, they expressed the operator explictly in a man- 
ner consistent with its boost weight decomposition, and then studied its eigen- 
bivector and eigenvalue structure in order to refine the alignment classification 
(i.e., to divide several of the algebraic types into subtypes). The hope of such a 
refinement is that it might be possible to prove more powerful results in one or 
more of these more restricted typesH 

From a computational point of view, a disadvantage of these two classifi- 
cation schemes is that one must solve a complicated set of equations to find 
the 'preferred' null frame (of vectors or bivectors) in which to do calculations. 
It would be very useful if there was a more constructive way of accessing the 
invariant classification information. In 4D, such an approach is given by the use 
of scalar curvature invariants (e.g., see [7]). In [2^, similar algebraic techniques 
were developed in higher dimensions using 'discriminants'. 

As an example, consider the bivector operator C defined in For a space- 
time of algebraically special type, the operator C will have a restricted eigen- 
vector structure [3T1 [TU]. That is, we can obtain necessary conditions for a 
particular algebraic type, in terms of the eigenvalue structure of the matrix. 

Studying eigenvalues of the bivector operator alone can only give so much 
information. However, as we will discuss in Section [2J there are many other 
operators that can be considered, constructed from various contractions of the 
Riemann tensor and its derivatives. It was shown in [20] that it is possible to 
derive necessary conditions for a spacetime to be of various algebraic types in 
terms of these other operators. 

Computing eigenvalues of large matrices explicitly is difficult, but in fact 
there is no need to do this in order to obtain the information that we require. 
In [2D] , two of the present authors showed how one can use an explicit algorithm 
to completely determine the eigenvalue structure of the curvature operator, up 
to degeneracies, in terms of a set of discriminants described in Section [2] and in 
the Appendix. Since the characteristic equation has coefficients which can be 
expressed in terms of the scalar polynomial curvature invariants of the curvature 
tensor, d these conditions (discriminants) can be expressed in terms of these 
polynomial curvature invariants. 

In particular, these techniques can be used to study the necessary condi- 
tions in arbitrary dimensions for the Weyl curvature operator (and hence the 

^In all of what follows, Greek indices a, /3, /i, i/, . . . are _D-dimensional spacetime indices. 

^ Note that the alignment types G and I are both of equal generality with respect to their 
possible (eigenbivector /eigenvalue) roots structure. It has been suggested for other reasons 
that the distinction between these types may not be significant in certain applications |19| . 

^ A scalar polynomial curvature invariant of order fc is a scalar obtained by contraction 
from a polynomial in the Riemann tensor and its covariant derivatives up to the order k. 
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higher dimensional Weyl tensor) to be of algebraic type II or D, and create 
syzygies which are necessary for the special algebraic type to be fulfilled. We 
are consequently able to determine the necessary conditions in terms of simple 
scalar polynomial curvature invariants in order for the higher dimensional Weyl 
tensors to be of type II or D. 

Let X be the set of all scalar polynomial curvature invariants. An important 
question is whether the X is unique for the spacetime under consideration. In 
particular, does there exists a continuous deformation, g^-, of the metric g so 
that Qr gives the same invariants as gl If such a deformation exists we will 
say that the metric is X-degenerate (otherwise, it will be said to be I-non- 
degenerate) . In 4D this question was discussed in [22] , and it was found that the 
only spacetimes that are not I-non-degenerate are the (very special) degenerate 
Kundt spacetimes. In particular, this implies that black hole spacetimes are X- 
non-degenerate and that the set of scalar invariants can be used to characterise 
such spacetimes. The same is believed to be the case in higher dimensions and 
it is I-non-degenerate spacetimes that will be the focus of study here. 

The purpose of the current paper is twofold. First, we demonstrate that the 
discriminant approach is reasonably practical in that even in the case of some 
very complicated metrics it is possible to compute the relevant discriminants and 
extract useful information from them, at least with the use of computer algebra. 
Unfortunately, for some of the metrics that we discuss we will find that doing 
computations for generic values of various parameters is too difficult, and hence 
we will be forced to study particular representative examples within families of 
spacetimes (for example, in the case of the Pomeransky-Sen'kov doubly spinning 
black ring [3]). 

Second, we look to better understand the links between these discriminants 
and the alignment type of spacetimes. Previous work [20] gave a set of necessary 
(but not sufficient) conditions on the discriminants for a spacetime to be type II 
(or more special). Here we illustrate these conditions in a number of examples, 
and also prove a partial result fProposition l2.5p that demonstrates that stronger 
conditions hold in the case of spacetimes with considerable isotropy. Given 
these additional conditions in the case of spacetimes with isometrics (at least 
for some set of points in the spacetime), we are able to see that in all cases 
discussed in this paper the necessary conditions relevant to the isometry group 
of the spacetime also turn out to be sufficient. This naturally leads one to 
speculate that there might exist a more complete result linking scalar invariants 
and alignment classification more closely. 

A useful application of these techniques to the alignment classification of 
spacetimes is also given. It is known [TB] that the singly-spinning black ring 
[3] is of type II on the horizon, but only of type I or G outside. The rele- 
vant discriminants were computed in [2D] and found to be consistent with this 
result. Therefore, it might be expected that the doubly-spinning black ring 
would also have this property. Here, we verify this by showing that the rele- 
vant discriminants are non- vanishing outside the horizon, [f] and hence that the 
doubly-spinning black ring must be of type G or I in the exterior region (it is 
easy to check, for example by explicit calculation in Gaussian null coordinates 
|23| . that all black holes are of type II or more special on the horizon). 

*A syzygy is an algebraic relationship between scalar polynomial curvature invariants. 
^To be precise, the discriminants are computed at a number of particular parameter values 
and for fixed values of one of the coordinates; see Section 5. 
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2 Review of the Discriminant Approach 

In this section we review the scalar invariant approach to the algebraic classifica- 
tion of spacetimes introduced in |20) , which links scalar invariants and bivector 
operators [TU] with the alignment classification [S]. The approach described in 
PO] was based on a set of algorithms for determining the structure of roots of 
a polynomial, derived in j24j [25j [26] (more details of these results are given in 
the Appendix). 

We will develop this work further in the case of spacetimes with isometrics, 
proving a new result that extends the usefulness of scalar invariant methods in 
this case. 

2.1 Curvature operators 

Given a curvature tensor, an operator describing an automorphism of some 
finite-dimensional vector space can be defined. The most well-known example 
of such a map is the bivector map ([3]) obtained from the Weyl tensor, mapping 
the n = D{D — 1 ) /2-dimensional space of 2-forms to itself. 

If the Weyl tensor is of a particular algebraic type in the alignment classifica- 
tion, then the associated operator C will have a restricted eigenvector structure 
[10] . One can extend this technique to include other operators, constructed ei- 
ther from the Weyl tensor, or from other curvature operators R G 7?., where TZ 
denotes the set of all curvature operators constructed from the curvature tensors 
and their polynomial invariants. For example, a second operator that can be 
constructed from the Weyl tensor is given by 

C:Y^,^ iC//Fp.. (4) 

acting on the D^-dimensional vector space of 2-tcnsors Y , while the tensor 

= C^^'-'^C^^.p (5) 

can be used to construct an operator 

T : X" T'^X'^ (6) 

acting on the _D-dimensional tangent space of the spacetime (in fact, it is often 
more useful to consider the operator constructed from the tracefree part of T) . 

2.2 Discriminant analysis 

For a given curvature operator, R, we can consider the eigenvalues of this opera- 
tor to obtain necessary conditions on various alignment types of the spacetime. 
In particular, requiring the algebraic type to be II or D will impose restrictions 
of the eigenvalues on the operator [20 . Consider the characteristic equation 

/(A) = det(R - Al) = 0. (7) 

This equation is a polynomial equation in A and the eigenvalues are the roots 
of this equation. The coefficients can be expressed in terms of the invariants 
of R. Therefore, we can give conditions on the eigenvalue structure expressed 
manifestly in terms of the polynomial curvature invariants (syzygies) of R. 
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The precise details of how this can be done are given in the Appendix. In 
the remainder of this section we briefly summarize the minimal information 
required to give a basic understanding of the language and notation used in the 
remainder of the paper. 

Given a curvature operator, R, acting on an n-dimensional vector space, we 
can construct a unique discriminant sequence for the characteristic polynomial 
/, consisting of n scalar functions on spacetime, that we will denote by 

[J^i,'fp2,...,I^„]. (8) 

We will make particularly regular use of the discriminant sequence for C in what 
follows, and for convenience we will abbreviate these discriminants by writing 
Di = ^Di throughout this paper. [f| 

Given a discriminant sequence (j66p . its sign list S is defined to be 

S = [signQDi), sign(^i?2), . . . , sign('f^i?„)] (9) 

(where sign(a;) = 1, 0,-1, respectively, as x > 0, a; = 0, a; < 0). 

Given any sign list S = [si,...,s„], we construct a revised sign list S as 
follows. First, we look for "internal zeros" of S (that is, for any subsequences 
of the form [si, 0, 0, . . . , 0, Sj], where Si ^ and sj ^0). If there are none, then 
we take S = S. Otherwise, we replace any such subsequence with 

S^, Si^ S^, Si, Si, Si, Si, Si, Si ... , Sj] (1^) 

in S. The revised sign list will therefore contain no "internal" zeros, but may 
have zeros at the end. 

This is useful because of the following result: 

Lemma 2.1 f |241 1251 120| ). Consider the revised sign list S of the discriminant 
sequence [^i, • ■ • , r^^ti]- Let K denote the number of sign changes and L 
denote the number of non-zero members of S. Then, the number of distinct 
pairs of complex conjugate roots is K , and the number of distinct real roots is 
L~2K. 

Hence, by applying the algorithm given in the Appendix to compute S (which 
can conveniently be done in an automated way using a computer algebra system 
such as, for example, Maple), we can learn a significant amount about the 
eigenvalue structure of R. 

2.3 Necessary conditions for Weyl type II/D 

As discussed in the introduction, for the Weyl tensor to be type II (or more 
special) the eigenvalues of the corresponding bivector operator C need to be of 
a special form. 

In fact, since the invariants of a type II spacetime are the same as for type 
D, we will assume for simplicity that the spacetime is type D for the purposes 
of the present discussion. In [10], it was shown that in this case we can choose 
a basis 

iArUi, iAn, rriiAmj, n A rrij (11) 



®Ref. 1201 used CHV to represent Dio, -Dg, -Dg; to avoid introducing yet more notation we 
do not do this here. 
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adapted to the vector basis ([1]). Hence, the space of bivectors can be split into 
three vector subspaces of boost weights +1, and —1, of dimensions {D — 2), 
{n — D + 2) and {D — 2), respectively. With respect to this basis, the bivector 
operator takes the form of an n x rt matrix {n = D{D — l)/2) 

C = blockdiag(Af, M*), (12) 

where M is a (D - 2) x (D - 2) matrix and \1/ is an (n - Z? + 2) x (n- D + 2) 
matrix |20|. Since the eigenvalues of M and M* are the same, it must be the 
case that: [j 

Lemma 2.2 ( |101 120] ). The Weyl bivector operator for a spacetime that is of 
type II or more special has at least {D — 2) eigenvalues of multiplicity (at least) 
2. 

Consider now the particular case of five dimensions. Here, the bivector 
operator acts on a vector space that is n = 10-dimensional and if the spacetime 
is type II then the bivector operator has 3 eigenvalues of (at least) multiplicity 
2. Hence the eigenvalue (Segre) type of the matrix is {(11)(11)(11)1111} (or 
more special, e.g. {(1111)(11)1111}). From the point of view of the discrimant 
analysis, it was shown in 20! that eigenvalue types consistent with at least 3 
matching pairs of eigenvalues can only occur if the last three discriminants in 
the sequence ([5]) vanish (see the Appendix for further details). Hence, we arrive 
at the most important result for this paper: 

Proposition 2.3 f[20j). If a five- dimensional spacetime is of alignment type II 
or more special, then 

Ds^O, Dg^O, Z?io-0. (13) 

These three equations are syzygies of order 90, 72 and 56, respectively. In 
principle, these can be computed in general, and written out in full in terms 
of the Weyl tensor. In reality, this is not practical, and it is far more useful to 
apply the algorithm used to construct the discriminants in the case of particular 
metrics. 

In this paper we shall also consider the operator T defined in ([5]) above. 
Following a similar approach to the bivector case, it can be shown that: 

Proposition 2.4 f[20|). For a five dimensional spacetime of algebraic type II 
or more special, the discrminants associated with the operator T satisfy: 

|7^5 = 0, fi^4>0, \D3>0, |i^2>0. (14) 

Note that jDr, = is a 40th order syzygy in the Weyl tensor (a 20th order 
syzygy in the square of the Weyl tensor), and hence is likely to be easier to 
calculate explicitly than the syzygies resulting from the bivector operator. 

^We have discussed the type D case here. In the type II case the block-diagonal form 
above reduces to a 'block lower-triagular' form for which a similar argument applies. The case 
of a spacetime that is more special than this (i.e., primary type III, N or O) is far simpler; 
this occurs if and only if the bivector operator is nilpotcnt. 
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2.4 Spacetimes with isometries 

Note that Proposition [53] gives necessary conditions for a spacetime to be type 
II or more special. These conditions are not sufficient. Indeed, the necessary 
conditions of Proposition 2.3 may be fulfilled in spite of the fact that a spacetime 
is of type G or I. Many examples of the non-sufficiency of this condition can be 
found by considering spacetimes with large amounts of symmetry; we will see 
several examples of this in the remainder of the paper. However, in this section, 
we will show how for spacetimes (or sets of points in spacetimes) admitting an 
SO{d— 2) isometry, there is a stronger necessary condition on the discriminants 
which may be of more use. 

More explicitly, it was shown in [TD] that a 5D spacetime with a Weyl tensor 
that has an 5*0(2) isotropy with spacelike orbits must automatically have 3 pairs 
of matching eigenvalues. Intuitively, this is because one can choose a basis in 
which the isometry makes all quantities invariant under SO{2) transformations 
mixing rn2 and 777.3, and hence both the 3x3 matrix M, and the 4x4 matrix 
must have two of their eigenvalues the same. Matching this with the result 
above, we see that one eigenvalue must be repeated at 4 times (since M has 
a repeated eigenvalue, and M* has the same repeated eigenvalue), and hence 
the eigenvalue type of the bivector matrix must be {(1111)(11)(11)11} or more 
special. This is now a basis independent condition. 

This condition can be interpreted in terms of the discriminants Di ; thus 

Proposition 2.5. // a 5D spacetime possesses a spatial 50(2) isotropy, and is 
of alignment type II or more special, then 

Dr^Ds^Dg^ Dio - 0. (15) 

Further details of this are given in the Appendix (of course the second half 
of the equalities follow trivially from Proposition 12.31 but the important point 
is that here they also hold as a result of the isometry itself, and hence we always 
need some additional conditions). 

Note that this is a pointwise result, and hence it also holds on lines of 
enhanced symmetry within a spacetime. The most important example of this 
is on an axis of rotation in an axisymmetric spacetime; we shall see several 
examples of this in the sections that follow. 

2.5 A caveat: Discrete symmetries 

Leaving aside the links to alignment classification for the moment, the bivec- 
tor methods [TU] also provide an independent algebraic classification in their 
own right; with algebraic types labelled by the eigenvalue (Segre) types of the 
bivector operator. 

Using the invariants only, we can only determine the eigenvalue type of the 
operator. For example, we cannot distinguish between the eigenvalue types 
{(1, 1)11...!} and {1, l...l(ll)}. This is a result of the fact that the spacetime 
is Lorentzian and that time and space directions have a different interpretation 
in spacetime. The discriminants cannot a priori distinguish the time and space 
directions. This is a caveat that is distinct to the Lorentzian nature of spacetime 
and is not present in analyses of Riemannian (i.e., signature (+++..+) ) spaces. 
Note also that this degeneracy is a discrete degeneracy, unlike the notion of 
I-degenerate metrics |22j which requires a continuous deformation. 
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3 The Sorkin-Gross-Perry soliton 

In order to illustrate some of these caveats, we now consider the example of the 
Sorkin-Gross-Perry (SGP) soliton, a class of 5-dimensional, vacuum solutions 
with a large amount of symmetry [27] . 
The 5D SGP metric is given by: 

\l + m/Rj \l + m/Rj 

where a = /3'^ + (3 + 1 (for a vacuum soliton solution) and we assume that 
m 7^ (so the spacetime is not flat; in the calculations below we set m = 1, 
which simply re-scales R), and the coordinates range over — oo < t, z < oo and 
m < R < oo. 

Generically, the solution does not contain a black hole; it closes off smoothly 
at i? = m. However, there are a number of special cases: 

(i) li (3 — 0, it is the Schwarzschild black string, which is of type D every- 
where. 

(ii) If /3 = 1, it is a boost invariant singular spacetime of type G. 

(iii) The limit /3 oo corresponds to the static Kaluza-Klein bubble; that is, 
the product of a flat time direction with the Euclidean 4D Schwarzschild 
solution, which is of type G everywhere. 

In other cases, these solutions are of Weyl alignment type G or I [55J [T5] . 
In fact, it was noted in p,9' that these provide examples of connected analytic 
spacetimes that are type G in some open region and type I in some other open 
subset, determined by the sign of a particular inequality. For example, if we 
choose parameters to > and /? = 1/2 in the above metric (which exhibits 
behaviour typical for the case < /? < 1 and m > 0) , then it is of type G for 
R ^ m but of type I for i? ^ m. Conversely, for (finite) /? > 1 and m > 0, it is 
of type I near R = m and of type G for i? ^ m) . d 

This spacetime is static: d/dt is a hypersurface orthogonal timelike Killing 
vector field. It was shown in [M] that all static spacetimes (and a large class of 
stationary spacetimes) in dimensions _D > 4 are necessarily of alignment type 
G, Ii, D or O. This implies that this spacetime cannot be of type II, III or N. 
Essentially, the reason for this is that static spacetimes allow for the discrete 
transformation 1 1— > — i which interchanges the two null directions. 

*By solving a quartic equation it was claimed in 28 that a null coframe exists in which 
the positive +2 (and negative —2) boost weight Weyl components vanish, and hence that the 
SGP metric is of type I everywhere. However, it was not noticed that the solutions for the 
roots of this quartic are not always real, and in some regions in which the roots are complex 
the algebraic type is in fact G (i.e., the preferred null vector is actually complex). Note that 
within the bivector formalism |10| these solutions are of equivalent algebraic specializations in 
the different regions. This is consistent with the results of 1191 discussed above and corrects 
the omission in |28l . 
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3.1 Discriminant Analysis 

We now look to compute the discriminant sequence for the bivector operator 
C for this metric. For all of the SGP spacetimes (with R > m), we find that 
Di = 10 (which merely sets a normalization) but also that 

Ds^Dg^ Dio = 0. (17) 

Hence, the necessary condition from Proposition 12.31 for the spacetime to be 
of alignment type II or more special is always met. However, this spacetime 
has sufhcient symmetry that the necessary conditions of interest here are the 
stronger ones of Proposition [53] In the general case, all other Di are non-zero, 
and these stronger conditions are not satisfied. 

To see this in more detail, consider a typical case P — 1/2 (and, without loss 
of generality, set m = 1). For convenience, define Q{R) = —7R^ — 7 + 6y/7R, 
which only has one real root with i? > 1 (namely Rq = (3a/7 + yT4)/7). We 
then find that, for R close to Rq, we have that 

D2 > 0, D3> 0, (18) 

and 

D4 = {+ve)Q\ Z?5 = (+we)Q^ = (+«e)Q^ D7 = {+ve)Q^°. (19) 

Hence, for R ^ Rq (i.e., Q 7^ 0) the sign list (and revised sign list) of 
this operator is [1, 1, 1, 1, 1, 1, 1, 0, 0, 0]. This contains no changes of sign, and 
hence Lemma |2. II tells us that the operator has no complex eigenvalues, and 7 
distinct real eigenvalues. Therefore, the eigenvalue type of the 5D Weyl tensor 
is {(11)(11)(11)1111}, corresponding to type I/G (consistent with the results of 
[19| and recalling that types G and I are of equivalent algebraic specialization 
within the bivector formalism [10 ). Note that even with D% — Dg = Diq — 0, 
this eigenvalue structure implies type G or I because the SGP spacetime has a 
spatial isotropy (as discussed above and in |21)). 

Note that on R = Rq, D4 — D5 — Dq = D7 = 0, and so on this line the 
spacetime could be of type D. This is indeed the case P^ . 

Let us next consider the special cases (i) — (iii). 

(i) In the case of the black string {(3 — 0), we have that 

Di = 10,^2 > 0,D3 > 0, (20) 

and 

Di = D5 = DG = D7 = Ds^DQ = Dw = 0. (21) 

This satisfies the necessary condition for a spacetime with an 5*0(2) isom- 
etry to be algebraically special and, in fact, the black string spacetime is 
of type D everywhere. 

ii) For /3 = 1 we have that Di > 0, for i = 1 — 5, and Di = 0, for i = 6 — 10. 
As noted above, this case is of type G. 

(iii) The static KK bubble spacetime /3 — > 00 also has 



Di = 10,L»2 >0,D3> 0, 



(22) 
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and 

Di^D5^DQ = Dj = D^ = D^ = Dio = 0, (23) 

and is of type G everywhere. 

For the black string and the static KK bubble, Di, D2, D3 are positive, while 
all other Di — 0, but the black string is of type D and the static KK bubble 
is of type G. However, there is a discrete symmetry (interchanging r and t) 
relating these two 5D spacetimes (which are related by a complex continuation) 
and hence their invariants are identical but their Weyl types are different (this 
is similar to the discrete, conformally flat, example given in [12] )• 



4 Supersymmetric Black Ring 

The results of this paper do not only apply to vacuum spacetimes. To emphasize 
this, we now move on to consider an interesting class of non-vacuum metric, 
corresponding to supersymmetric, rotating black holes in five dimensions. The 
simplest example of such a spacetime is the extremal, charged, rotating BMPV 
black hole |29| in 5D minimal supergravity, which has a horizon of spherical 
topology. It is already known that the BMPV metric is generally of Weyl type 

Later, a supersymmetric black ring solution of the same supergravity theory 
was presented in [30], with a black hole horizon of topology x S^. In fact, 
this solution can be seen as a generalization of the BMPV metric, which can be 
obtained as a particular limit of the supersymmetric black ring family. 



4.1 Metric and properties 



Since the BMPV metric is a particular case of the supersymmetric black ring 
(having a horizon topology x S'^), one would expect that the supersymmetric 
black ring is not of Weyl type II or more special, but this has not been proved. 
This question is something that we can attack with relative ease using these 
discriminant methods. 

The line element of the supersymmetric black ring is 



ds^ = -fix, yfidt + ujf + fix, y)-^ds\n^), 



where 



= 1 



ix-y)~^ix^-y^) 



U! 



uj^d(j) + oj^dtp, 

-s^il-x^)[5Q-q^i3 + x + y)], 
§9(1 + y) + 8^2(1 -2/')[3Q- 9^(3- 



and the four dimensional flat space metric is written as 



ds2(R4) 



(x - yy 



dy' 



1 



+ iy'^ ~ l)d^P^ + 



dx^ 



1 



+ (1 - x^)d(l)^ 



(24) 

(25) 
(26) 
(27) 
(28) 

(29) 



Admissible coordinates values are — 1 < x < 1, ^00 < y < ^1 and (jjjip are 2tt- 
periodic; and the parameters q and Q satisfy q > and Q > (which implies 
that /(x,2/) >0). 
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The black ring horizon hes at y = — oo in these coordinates. Of course a 
different choice of coordinates can be made that includes the black hole horizon; 
and typically one must do this when trying to do any computation of a quantity 
on the horizon. A great advantage of scalar invariant approaches is that this 
is not necessary; one can compute a scalar invariant everywhere outside the 
horizon and then extend it onto the horizon itself by continuity. Hence, we are 
able to work in this relatively simple coordinate system throughout. 



4.2 Discriminant Analysis 

Unfortunately, computing the discriminants Di explicitly for arbitrary values of 
the parameters has proven to be difficult; hence we are forced to pick particular 
values of the parameters (i.e., particular spacetimes in the family) to analyse. 
The parameter R merely fixes the length-scale, and hence we can set R = 1 
without loss of generality. We have then computed the discriminants for various 
values of the parameters q and Q; for the purposes of the present discussion, we 
consider the case q = 1/2, Q = 9/4. [f| 

Let us first consider the Weyl operator on the horizon y — ~oo. For the 
Weyl bivector operator C, the discriminants can be evaluated, and give 

Di = 10, L»2 = 180, Ds = 4608; D4 = D5 = Dq = D7 = Ds = = fio = 0. 

(30) 

Hence, this is consistent with the black ring being type II or D on the horizon, 
as we know it must be. Explicitly, the eigenvalue type of the Weyl bivector 
operator is of the form {1(111)(111111)} on the horizon, while the operator T 
is of eigenvalue type {(11)(111)} there. 

Outside the horizon, the computations are rather more complicated. Com- 
puting the operator T indicates that the discriminant jD^ takes the form: 



_ l)2(y + 1) 2(^_1)2(^ + 1)2^(^^ y)154 

{x - yyo 



\D,^F{x,yy-^ '-^ ' \^ '\ ' , (31) 



where F{x,y) is some (known) polynomial in x and y and / is the metric function 
defined in ([25)1 . Recall that wherever jD^ is non-zero, the Weyl tensor must be 
of type I or G. 

Checking that F{x, y) is non-zero is relatively difficult. We have verified that 
this is true for the following cases (among others): 

• X ^ 0, y general. 

9 X — ±1/2, y general. 

• y ~ ~2, X general. 

This shows that the metric is of type I/G except possibly for a few special 
values of {x,y). For the SUSY black ring, some of the special cases are as 
follows: X — y, which represents a curvature singularity, occurs only at x=y=-l 
(for the coordinate ranges used here) and corresponds to asymptotic infinity. 
X = ±1 is the plane of the ring (the axis of (j) rotation) and y = — 1 is the axis 
of T/j-rotation. 



^ We also computed the discriminants in the cases q = 1/4, Q = 1; q = 1/2, Q = 1/2; 
= 1, Q = 2. 
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5 Doubly-spinning black ring 



Pomeransky & Sen'kov [1] constructed a family of exact solution to the vacuum 
Einstein equations in 5D, corresponding to a black ring with both angular mo- 
menta non- vanishing. The algebraic type of this metric is currently unknown. 
Since this family contains the Emparan-Reall black ring l3l as a special case, 
it seems unlikely that other doubly-spinning black rings are more algebraically 
special than this solution. To prove this directly, one must demonstrate the 
non-existence of solutions to the appropriate alignment equations; the difficulty 
in doing this lies in the complexity of the metric. 

Therefore, as a demonstration of the utility of the discriminant techniques, 
let us do this by another method, and using the discriminant techniques dis- 
cussed above prove that the doubly spinning black ring is not type II or more 
special, except on the horizon. Unfortunately, we are still not able to do this 
explicitly for general parameters R, A, (defined below), as the resulting equa- 
tions are too complicated even with computer algebra, but we are able to do it 
for a representative sample of different parameters covering all known cases of 
qualitatively different behaviour. 



5.1 Metric and parameters 

The metric can be written in the form |31) : 



H{x,y) 



{dt + 0)2 



dx^ 



dy' 



G(x) G(y) 



+ 



A{y, x)d(tP' - 2L{x, y)d(l)d-)p - A{x, y)d'ip'^ 
H{x,y)H{y,x) 



(32) 



where the coordinates lie in the ranges — oo < t < oo, < (p,ip < 27r, —1 < x < 
1,— oo<y<— 1. The metric is globally asymptotically flat, although this is 
not manifest in these coordinates, where asymptotic infinity corresponds to a 
point {x,y) = (-1, -1). 

The one- form characterizing the rotation is given by = rt^{x,y)d'ip + 
fl^{x,y)d(j), where 



0?/ 



i?A^2((l + I/)2 - A2) 1 + y 



H{y,x) 



1 - A + i^ 



(l + A - J/ + x^yv{l - A - I/) 



2vx{l~y)) (33) 



and 



(34) 



14 



A. Coley, S. Hervik, M. Durkee and M. Godazgar 



while the polynomial functions G, iJ, A, L read 

G{x) = {l-x^){l + Xx + iyx^), (35) 
H{x, y) = l + -v"^ + 2\v{l - x^)y + 2xX{l - y'^v'^) 

= Wy'^v{l-\^ -v"^), (36) 
L{x,y) = \^/D{x - y){l - x^){\ - y'^)[l + - v"^ + 2{x + y)\v 

-xyu{l - \^ - y-")], (37) 
A{x,y) = G{x){l-y^)[[[l-yf -X'){l + u) + y\{l-\^ + 2u--iu^)\ 
+G(y) [2A2 + a;A((l - vf + A^) + ^^((1 - vf - + u) 

+x^\{l - A^ - iv^ + 2v^) + x'^v{l - v){l - A^ - v^)] .(38) 

The (outer) horizon lies at G{y) = or 



y^y-^ — 2. — ■ ^^^^ 

In this paper, we will not consider the black hole interior. 

The parameter i? > sets the length-scale in the spacetime, while the di- 
mensionless parameters A and v are restricted to 

i? > 0, < < 1 , 2^/D <X<l + v. (40) 

When X = V = Q the metric reduces to that of flat spacetime, written in ring- 
like coordinates. The singly spinning limit is obtained by setting v ~ 00 The 
bound A > 2^Jv corresponds to a Kerr-like upper bound on the rotation of the 
black ring around the S^. When A = 2y/v, the horizon is degenerate and we 
obtain an extremal black ring. 

It was shown in [33l [31] that the properties of the black ring (outer) ergo- 
surface vary drastically with the values of A and v. For 2^/1) < A < I — 
the ergoregion has the same topology as the event horizon (S*^ x 5^), but when 
I — u<\<l + v the ergosurface has topology 5^ U 5^, i.e. it consists of the 
union of two disjoint spheres, one surrounding the entire horizon, and the other 
lying in the centre of the ring. There is a critical case at A = 1 — where the 
ringlike ergosurface just 'pinches off' on an . 

There is a generalization of this family that allows for a 'non-equilibrium' 
black ring, supported by a conical singularity. In this case, the metric functions 
take an extremely complicated form in the general case in which the black ring 
is not in equilibrium [34j . The extremal MP solution can be recovered as a limit 
of the extremal, non-eqilibrium solutions in this family, with — > 1, A — > 2. 



5.2 Discriminant Analysis 

We now look to evaluate the discriminants Di and jDi for this metric. Unfor- 
tunately, computing the determinants for general values of the parameters A, 
V, and the coordinates x and y, turns out to be too complicated, even with the 
use of Maple. 

^"^To obtain the precise form of the singly spinning black ring given in the review articles 
[T]|32], one must identify B? = 2A;^(1 -|- A)^ and rename A — > 1^. 
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We will therefore compute them for particular solutions within the black 
ring family. There are (at least) five distinct parts of parameter space (A, v) 
that give some kind of known, qualitatively different behaviour for the solution. 
To try to ensure that we do not miss any important behaviour, we include an 
analysis of each of these cases: 

• Case 1. 0<A<l,j/ = 0: Singly spinning, non-extremal, x 5^ 
ergosurface (e.g. A = 7/9, v = 0). 

• Case 2. 2^/// <A<1 — i^, i/>0: Doubly spinning, non-extremal, x S"^ 

ergosurface (e.g. A = 7/9, u = 1/9). 

• Case 3. \ = 1 — > 0: Doubly spinning, non-extremal, 'pinched' 

ergosurface (e.g. A = 8/9, = 1/9). 

• Case 4. 1 — j/<A<l + i/, i/>0: Doubly spinning, non-extremal, 

U ergosurface (e.g. X = l, i/ = 1/9). 

• Case 5. A = 2y/u, v > 0: Doubly spinning, extremal (e.g. A = 2/3, 
u=l/9). 

First, we consider the discriminants Di. Even having fixed A and v, it is not 
possible to evaluate the discriminants for general x and y. Hence, we will first 
fix X (separately choosing x = and x = 1/2), and vary y in order to capture 
the change of behaviour both on the horizon y = yn and on the '0-axis y = —1. 

Doing this for both values of a; = and x = 1/2 , we find that for the first 
four cases above: 

£>8 cx (2/ + lf{y - ynfiy - xY^\ <x {y + l)\y - yh)\y - xf^\ 

D^oo,{y + lf{y-yhf{y-xY^\ (41) 

where only the relevant factors have been kept. That is, we have only included 
factors that appear in all the three discriminants above. Recall that all three 
discriminants must necessarily vanish for the solution to be type II at that 
point. We have also neglected factors that vanish for values of y outside its 
range. 

For the final (extremal) case we find for both values of x = and a; = 1/2 
that: 

D^<x{y + lf{y - yh)'%y - a:)"^ D, <^ {y + l)\y - y^f^y - xf^\ 

D^^o,{y + lf{y~yhr{y-xf^\ (42) 

and also 

D^^{y-ynf{y-xf\ o, {y - ynf {y - xf\ (43) 

Now, wc fix y depending on the value of A and v such that it is outside the 
horizon. For the first case where the black ring is singly spinning, the ergosurface 
is at a fixed value of y. For our example, A = 8/9, u = 1/9, y^ = —65/63. 
Thus, we choose three values of y inside the ergoregion (y = —6/5), outside the 
ergoregion {y = —65/64) and on the ergosurface {ye = —65/63) and calculate 
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the discriminants to capture their x— dependence. We find that for all three 
values of y 

Dso^ix- Ifix + l)Hy ~ xY^\ D.ocix- l)i°(x + l)i°(y - x^^^ 

i^iooc (x-l)i«(a; + l)i6(y-:r)i«", (44) 

and also 

De(xix~lfix + lfiy-xf", Dt (x {x - l)\x + l)^ {y - xf\ (45) 

For the next four cases for which the black ring is doubly spinning, the 
ergosurface is not at a fixed value of y. Thus, in each of these cases we pick two 
values of y, one inside the ergoregion and the other outside: 



• Case 2 (A 


7/9, 1^ 


- 1/9): 


V 


^ -21/20, y^ 


-3/2. 


• Case 3 (A = 


8/9, 1^ 


= 1/9): 


y 


= -41/40, y = 


-6/5. 


• Case 4 (A = 


l,v = 


1/9): y 




-161/160, y = 


-11/10 


• Case 5 (A = 


2/3,1^ 


= 1/9): 


y 


= -11/10, y = 


-2. 



In all of the cases above, for both values of y, we find that 

DsOc{x- ifix + ifiy - :r)"2, i^g oc (x - l)^ix + ifiy - xf^\ 

Dw(x{x-lfix + lf{y-xy^". (46) 

In summary, the analysis of the discriminants Di suggests that the black 
ring can only be of type II or more special at 

y^yh, y = -l, x^±l; (47) 

i.e. at the horizon (where we already know it to be type II or more special 
[23 [H]), the ■(/'-axis {y = —1) and the 0-axis (a; = ±1). 

Finally, we consider the discriminants jDi. To make the analysis easier we 
use the results obtained from the analysis of the discriminants Di. Recall that 
the necessary conditions on jDi for the solution to be type II or more special 
is that: 

^D2>0, \D3>0, \Di>0, |i?5-0. (48) 

As expected, we find that these conditions are satisfied for y = yh for all five 
cases. However, for the axes {y — —1 and a; = ±1) these conditions are only sat- 
isfied for the singly spinning black ring. For the other cases, jD^ — 0. However, 
one of the other discriminants may be negative. Thus, the axes only satisfy the 
necessary conditions for type II or more special for the singly spinning black 
ring. 

Therefore, we conclude that the double black ring is not type II or more 
special except at the horizon. In addition, for the singly spinning black ring, 
the axes, y = —1 and x = ±1, also satisfy some of the necessary conditions. 
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6 Other examples 

There are many other higher-dimensional spacetimes of interest. As a further 
test of the usefulness of these techniques, we will briefly study what discriminant 
analysis gives us for a collection of further spacetimes, that look rather different 
to the asymptotically flat spacetimes considered so far. In doing so, we hope 
to build up further intuition for whether discriminant analysis is likely to be a 
useful tool for studying higher-dimensional spacetimes. 

6.1 Black holes with nilmanifold horizons 

A 2-parameter family of exact vacuum solutions representing a 5D black hole 
spacetime with a local 3-dimensional nilmanifold horizon was found in [36) . The 
metric can be written as 

+ ^8/3 (^dz - ^-xdy^ (49) 

Using the operator T, we can determine the type outside the horizon. The 
spacetime admits a spatial isotropy, which the operator T will inherit; and con- 
sequently jZ?5 = 0. We can simplify expressions by fixing r = 10 Computing 
jZ?4, we get 

cx (2211M + 104L2)2(-934241Af2 + 3102L2A/ + 312L'')2 (50) 
x(432575M3 - eilOSL^A/^ - 1320L^M -f UAL^f 
X (6320919M3 - lASAh^QL^M"^ + 30096L^M -f 288L^)2 
X (3955732M3 - 418055^^7^/2 - 27324L''M + iiQL^f 
X (-251801M2 + 3234L2Af + 234L'')2(-11M + 2L^f. (51) 

When M — 2^^/11, the horizon lies at r = 1, and hence \Di = 0. Away from 
the horizon this discriminant is generally non-zero (except for a few special 
values). The eigenvalue structure is thus {1, 11(11)}, and hence the Weyl tensor 
is of type I/G outside the horizon. 
On the horizon, we get 

\Dr, = \Di^\D^^Q, \D2^0, Fl = 0; (52) 

consequently, T has eigenvalue type {(1111)1} on the horizon. This is consistent 
with the spacetime being of type D on the horizon. 

^'^We can do this without loss of generality, as we can recover the r-dependence if we note 
how the scalars and metric scale with r. If I{L, M) is an nth order invariant, then we can 
recover the r-dependence as follows: 



I{L,M) ^ r"3"/(Lr2/3,Mr-^/3). 
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6.2 Black holes with solvmanifold horizons 

A black hole with a local 3-diniensional solvmanifold horizon [37] was also found 
in |36] . The metric can be written 

(53) 

The asymptotics of this spacetime is locally that of a homogeneous Einstein 
solvmanifold [551 ISH] ■ Such spaces are generalisations of AdS spacetimes (which 
is the simplest form of a solvmanifold) and are locally distinct from such. In- 
terestingly, these black hole spacetimes share many of the features of topologial 
black holes in AdS but are topologically different. 

This has a discrete symmetry {x,y,z) i— )■ {—x,z,y) which interchanges two 
of the spatial eigenvalues of the operator T. Consequently, yiDs — 0. However, 
computing the discriminant jD^ (setting r = 1) and requiring M > gives: 

\D4 = (+ve){L^ ~9M)^. (54) 

Requiring positive mass, i.e. M > 0, this is only zero at the horizon, = 9M; 
consequently, this metric is of type I/G outside of the horizon. 

On the horizon, this degenerates to type {(11)(111)}. The Weyl operator 
can easily be computed on the horizon and it is easily verified that it is of type 
D. 

6.3 5D Generalised Collinson-Prench spacetimes 

This is not a black hole metric but rather a one-parameter family of cosmological 
solutions in 5D [ID]. The family of solutions has the metric: 

ds^ = _di2 ^ ^2 ^^2 ^ ( t{l~pfs^~V6(l+p}r.^y ^ ^^^^^ ' 

+ i6(i+P)''e4^™dz2 + ffip(i+p)«e''P^"du;^ (55) 

where 

yi±Z±Z 1 (56) 

For p = 1, this has an isotropy, so we expect jD^ ~ here. To calculate 
the discriminants we note that this spacetime is self-similar with a time-like 
homothety. For an nth order invariant, the t-dependence will be t~'^"; therefore, 
there is no loss of generality to set t = 1 in the final expressions. In general 

(t = l): 

cx p2°(p-l)2(l+p)'^(l+p + p2)(57p^"-8l/ + 144/ + 94/ 
-121/ -I- 390p^ - 121p^ + 94/ + 144p2 - sip + 57) 
X (34/ + 17p^ + llp^ -ip + lf{7p^ - 3/ + 1739^ + 17p2 + 34)2. i^r^^^^ 



^^Here, if / is an nth order invariant, we can recover the r-dependence via 

I(L, M) ^ r-^"I{Lr, Mr'^). 
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Therefore, so long as p is different from 0, ±1 (and some other very special 
values), the spacetime is of type I/G. 

For p = 0, for example, jD^ = 0. Indeed, for p — all invariants are zero; 
hence, it is a VSIq space [9]. |ij Since the case p = is VSIq, we can consider 
the derivative operator: 

r~ a /^bdef:a/^ ^ a /^bdef:h/^ ^t^Q\ 

c — ^ ^bdef ;c — c'-' t'bdef;h- (OO) 

Computing the discriminants in the case p — 0, we get: 

|D5 = \Di = 0, pi > 0, and = 0, (59) 

which shows that the eigenvalue type is {(111)11}. Therefore, since this metric 
is VSIq but not VSIi , this metric is I-non-degenerate [22 • 

7 Discussion 

The study of higher dimensional black holes is of current interest. The algebraic 
classification of spacetimes has played a crucial role in understanding black holes 
in 4D, and it is likely to play a similar role in higher dimensions. There are a 
number of algebraic classifications of spacetimes in higher dimensions (utilizing 
alignment theory, bivectors and discriminants). Previous work has produced 
a set of necessary (but not sufficient) conditions on the discriminants for a 
spacetime to be of a particular (special) algebraic type. We have demonstrated 
that this discriminant approach is a reasonably practical approach in that even 
in the case of some very complicated metrics it is possible to compute the 
relevant discriminants and extract useful information from them. Indeed, in 
this paper we have illustrated these techniques in a useful application to the 
alignment classification of a number of important spacetimes. In particular, we 
investigated the Sorkin-Gross-Perry soliton, the supersymmetric black ring, the 
doubly-spinning black ring, and some other higher dimensional spacetimes. 

For the Sorkin-Gross-Perry (SGP) soliton spacetimes all of D^, Dg, and Diq 
vanish. In the general case we have that all other Di are non-zero. Note that 
even with Dg = Dg — Diq — 0, this eigenvalue structure implies type G or I 
because the SGP spacetime has a spatial isotropy. These SGP examples thus 
illustrate the application of Proposition 2.5 in the discrimination analysis. Note 
that on R = Rq, D4 = — Dq — Dj = 0, and so on this line the spacetime is 
of type D. 

We next considered the 5D supersymmetric rotating black (SBR) holes (that 
include the extremal charged rotating BMPV black hole of [2£,). Since Ds — 
Dg = Dio = 0, this is a signal that it is of type II/D on the horizon. Explicitly, 
the eigenvalue type of the Weyl bivector operator is {1(111)(111111)}, while the 
operator Tj' = CcdeaC'^'^'^'' is of eigenvalue type: {(11)(111)}, on the horizon. 
As long as jD^, is non-zero, the Weyl tensor is of type I/G. 

As a demonstration of the utility of the discriminant techniques, we also 
proved that the doubly spinning black ring (DSBR) is not of type II or more 

A VSIi space is a spacetime for which all of the zeroth to i-th order curvature invari- 
ants (i.e., the scalar invariants constructed from the Ricmann tensor and its first i covariant 
derivatives) vanish [9]. 
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special, except on the horizon. Indeed, we first showed, using the discriminants 
Di, that the black ring can only be of type II or more special at the horizon 
and on the ■i/^-axis and the 0-axis. We then considered the discriminants ji'i, 
and showed that the necessary conditions for the solution to be of type II or 
more special on these axes are only satisfied for the singly spinning black ring. 
We concluded that the double black ring is only of type II or more special at 
the horizon. 

There are many other higher dimensional spacetimes of interest. We also 
analysed the black hole solution with a Sol-horizon and the black hole solution 
with a Nil-horizon and the 5D generalised CoUinson- French solution. 
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A Discriminant Analysis: The Algorithm 

At a point in spacetime, any curvature operator R is a linear map from one vector 
space to another. In this Appendix we describe an algorithm for determining 
the eigenvalue structure of this linear operator and, in particular, give criteria 
for different Segre types. This is a review of the algorithm described in |20j . 
based on the work of [Ml HSl US ; to analyse the characteristic polynomial of the 
operator. 

The criteria for various eigenvalue types will be given in terms of a set of 
'syzygies' (scalar polynomial invariants) which can be used to characterise the 
various eigenvalue cases; i.e., they are discriminants. For our purposes, for an 
operator R acting on a n-dimensional vector space, the relevant discriminants 
will be the ten spacetime scalars I < i < n. 



A.l Discriminant sequences and primary syzygies 

More explicitly, we can write the characteristic equation as a polynomial 

/(A) = det(Al - R) = aoA" + aiA^-^ + ... a,A""' + . . . + a„ = 0. (60) 

Choosing a normalisation such that oq = 1, the other coefficients can be 
expressed explictely in terms of the polynomial invariants 

i?i=Tr(R), i?2 = Tr(R2), i?3 ee Tr(R3), etc, (61) 
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of R, b^F 
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Ri 




^3 R2 


i? 



Next, define the (2n + 1) x (2n + 1) discrimination matrix 



(62) 
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On-l 

























nao 


(n — l)ai 











.0 














ai 




an. 






















(63) 


or more explicitely 
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(65) 










I 0, 






n + i + 2 < j 


<2n+l 







We are now in a position to make the main definition of this Appendix: 

Definition A.l. Let R be an operator with characteristic polynomial /, map- 
ping an n-dimensional vector space to itself, as above. Then, gZ?i is defined 
to be the determinant of the submatrix formed from the first 2i rows and 2i 
columns of Disc(/). 

The discriminant sequence of the polynomial / is the sequence 



(66) 



When rewritten in terms of the curvature invariants {i?i,i?2, • • •}, using (|62p. 
these discriminants are the primary syzygies for the operator R. 

Note that if — then there can be at most {n — 1) distinct eigenvalues, 
and hence there must be a repeated eigenvalue. More generally, if ^n-j = 
. . . = ^„ = 0, then there are at most n — j distinct eigenvalues. We can 
use this to find necessary conditions for the an operator to have a particular 



^*Note that if an = 0, then the eigenvalue equation trivially factorises and we have a zero 
eigenvalue. Furthermore, if a„ = a„_i = . . . = Qn— fc = Oi then there exists a zero eigenvalue 
of multiplicity fc + 1. Therefore, before beginning the algorithm one should check for the 
existence of zero eigenvalues, and work with the reduced polynomial f{X)/\^. 
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algebraic type in the alignment classification. In particular, this can be used to 
determine the number of real roots of the characteristic equation, as described 
in Section [2] 

Multiple factor sequence. For polynomials of order 3 or less, the discrim- 
inants are sufficient to determine the complete eigenvalue structure. How- 
ever, for higher-order polynomials (corresponding to operators acting on higher- 
dimensional vector spaces), we need more information. 
To get this, consider the following definition: 

Definition A. 2. Let M{k, I) be the submatrix of the discriminant matrix 
Disc{f) formed by the first 2k rows and the first {2k — 1) columns + {2k + /)th 
column and construct the polynomials: 

k 

Afc(/) =;^det[M(n-fc,i)]x^-\ (67) 

for fc = 0,l,...,n— 1. The multiple factor sequence of the polynomial / is 

[Ao(/),Ai(/),...,A„_i(/)] (68) 
This sequence is useful because 

Lemma A. 3 f |24[ I25j). // the number of zeros in the revised sign list of the 
discriminant sequence of f{x) is k, then the greatest common divisor of f and 
its derivative f is given by Afc(/) = g.c.d.{f{x),f'{x)). Hence, the g.c.d. of f 
and /' is thus always in the multiple factor sequence. 

Indeed, the polynomial A(/)(x) = g.c.d.(/(a;), f'{x)) is the repeated part of 
f{x), because if A(/) has k real roots of multiplicities ni, n2, ...,71.^, and / has 
m distinct real roots, then / has k real roots of multiplicities ni -|- 1, n2 -I- 1, 
-f 1, and m— k simple real roots (similarly for complex roots). Therefore, 
by considering A(/) we reduce the multiplicities of all the roots by 1. 

We can now consider the discriminants of the polynomial A(/) in the same 
way as we computed the discriminant sequence of /. We will call the dis- 
criminant sequence of A(/) {"i^i, "i?2, "^^3, "-Bfe}. We can now use these to 
determine the revised sign list of the ii'-sequence, etc. We can repeat this proce- 
dure and consider A(A(/)) = A^(/), A'^(/) etc (we label the next discriminant 
sequence F in the table below). 

Iterating this, we eventually arrive at some j for which the revised sign list 
of A^ (/) contains no zeros. We can then compute the number of real/complex 
distinct roots of A^ (/). This in turn allows us to determine the roots and 
multiplicities of A^~^{f), which again enables us to determine the roots and 
multiplicites of A^~^(/) etc. At the end of this process, we have a complete 
root classification for f{x). Further details of this procedure are described in 

Using these methods, it is possible to give necessary and sufficient conditions 
on the various discriminants for an operator to have a particular eigenvalue type. 
From these conditions, we can extract the necessary conditions on Dg, Dg, Diq 
etc that we make extensive use of in this paper. 

For the 10-dimensional trace-free Weyl operator, a partial table demonstrat- 
ing this is as follows: 



Classification of 5D spacetimes using scalar invariants 



23 



^10 


Dq 


^8 


Dj 






10^3 


Eigenvalue type 


^0 














{111..1} 





^0 












{(11)11-1} 








^0 










{(ll)(ll)l...l} 








^0 











{(111)1-1} 











^0 






^0 


{(ll)(ll)(ll)l..l} 











^0 




^0 





{(lll)(ll)l..l} 











^0 










{(1111)11..!} 
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